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ABSTRACT 

We report initial results from a Herschel program to search for far-infrared 
emission from cold dust around a statistically significant sample of young brown 
dwarfs. The first three objects in our survey are all detected at 70/im, and 
we report the first detection of a brown dwarf at 160/im. The flux densities are 
consistent with the presence of substantial amounts of cold dust in the outer disks 
around these objects. We modeled the SED's with two different radiative transfer 
codes. We find that a broad range of model parameters provides a reasonable 
fit to the SED's, but that the addition of our 70/im, and especially the 160/im 
detection enables strong lower limits to be placed on the disk masses since most 
of the mass is in the outer disk. We find likely disk masses in the range of a 
few xlO~^ to 10~^ Mq. Our models provide a good fit to the SED's and do not 
require dust settling. 
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Introduction 



Young brown dwarfs exhibit circum- "s t ellar" disk phenom ena much like the i r mor e 



massive counterparts, e.g. iKlein et al.l (|2003[ ). lApai et al.l ( l2004j ) and iLuhman et al.l ( I2OIOI ). 
Although there are quantitative differences in physical parameters for disks arouiid sub - 
stellar objects such as the gas chemistry and degree of dust processing ( iPascucci et al.ll2009l ). 
the distribut i on of properties is relatively continuous across the sub-stellar boundary, e.g. 
Scholz et al.l ( 120091 ). Brown dwarfs (BD's) provide a qualitatively and quantitatively dif- 
ferent physical environment in which to study disk structure and evolution, particularly 
in conditions more favorable to lower mass planet formation and with lower stellar heat- 
ing and ionizing fluxes. Observational selection effects, however, make the study of BD 
disks more difficult than for T Tauri stars because of the low luminosities. The quantity 
of cold, T < 150K, dust in the disks is difficult to determine without sensitive measure- 



ments at A >30/xm. Ther e are two reported Spitzer detections of BD's at 7 0um (|Guieu et al. 



200/ 



Riaz fc Gizis 2008), and severa l have been detected at A ~ 1 mm ( Klein et al 



2003 



Scholz. Jayawardhana fc Woodll2006l ) with modest S/N. But the detected number is painfully 
small for any statistical investigation of the cold dust mass, disk flaring, dust emissivity, grain 
growth, and other properties. 



We describe our earliest results from a Herschel (IPilbratt et al.ll2010|) GTl (Guaranteed 
Time, Phase 1) program with the PACS photometer (jPoglitsch et al.ll2010l ) that will eventu- 
ally provide sensitive photometry at 70 and 160/im of 50 young BD's in several star-forming 
regions spanning a range of ages and spectral types. Our earliest sample is a strong function 
of Herschel scheduling and comprises three objects, two in the relatively older (8-10 Myr) 
TW Hya association (TWA), and one in the younger (~ 2 Myr) Chamaeleon I region. 



2. Observations and Data Reduction 

All the observations used the PACS "mini-scan-map" mode with similar integration 
times. This mode simultaneously images at 70 and 160/im with two scan maps at the 
recommended relative angle for good 1/f noise reduction and high sensitivity over an area 
~ 60"x90". Table [U lists the relevant parameters and Table |5] lists the AOR's. 

The data were first processed with the Herschel Interactive Processing Environment 
version 7.0, with standard high-pass filtering for point-source observations to produce fits files 
of the image, coverage, and uncertainty. The uncertainty images are not yet reliable, so we 
estimated uncertainties as described below. Figure [T] shows the final mosaicked images. We 
also processed each AOR separately to test that the individual observations were consistent 
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with the combined set for each object, i.e. the expected a/2 reduction in S/N. 



The final proces sing utihzed the psf-fitting photome try tool, c2dphot, from the c2d 
Spitz er Legacy Team (iHarvev et al.ll2006l : lEvans et al and based on the earlier DOPHOT 

tool ( Schechter. Mateo, fc Saha 19931 ). This tool can be used in various modes, including 
finding peaks above the background and fitting a psf to the local maxima, or fitting a psf 
to a fixed position, a mode useful for estimating noise and determining upper limits. Our 
quoted noise was determined this way by fitting a half-dozen arbitrarily chosen, nominally 
empty parts of the image within the high-coverage area. Upper limits for the two objects 
undetected at 160/im were also determined this way, and also by inserting artificial sources in 
the images at various fiux levels. This tool produces aperture fiuxes for comparison with the 
psf-fitted values, and those values agree well. In addition, aperture-fiux "curves-of-growth" 
show good agreement with the psf-fit photometry. 



Observational Results 



Figured] shows all three objects clearly detected at 70/im and SSSPM1102 also at 160yum. 
The level of structured background at 160yum is higher for the other two BD's; the upper 
limits at 160/im for those two were set by this structure rather than by the instrumental 
sensitivity and DC sky background. There is some weak, diffuse emission at 160/im close to 
the position of 2M1207, but the offset from both the nominal source position and the 70/im 
centroid is probably too large to be consistent with emission from the object. In addition to 
our confidence in these detections based on inspection of the images, we note that c2dphot 
reliably extracted the sources as bona fide point sources with no confusion in its most general 
"source-finding" mode. The small offset from the nominal position for SSSPM1102, ~ 2", 
is essentially identical in both wavelength channels and well within the typical Herschel 
pointing uncertainty. 



The derived flux densities and upper limits are 



also been observed by Spitz er with IRAC and MIPS (IRiaz. Gizis fc Hmielll2006l : iRiaz fc Gizis 



spectral range (IRiaz fc Gizis 



isted in Table 121 All three objects have 



2008: Luhman et al. 2010). Two were observed with the IRS instrument over its whole 



2007 



2008 



Morrow et al.ll2008f l. and one (IS0138) over the 



lO/xm sihcate feature (jPascucci et al.ll2009l ). All are detected in the 2MASS catalog and 
have reasonably well-determined spectral types as described below. We have collected these 
and available shorter wavelength data into the SED's in Figure [2l 

Figure [2] shows that the SED's of the two BD's in the TW Hya association are much 
more similar to each other than to IS0138 in Cham I. In particular, the [8]-[24] color and 
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to a lesser extent the [8]-[70] color of IS0138 is substantially redder than for the other two. 
The two TWA BD's also show strong excess emission to wavelengths as short as the IRAC 
3.6/xm band. Most importantly, all clearly show emission in the far-infrared indicative of 
substantial amounts of cold dust. 



4. Modeling 



We have modeled the SED's of these three BD's with two different codes, MC3D 



described by IWolf Herining fc StecklumI (119991 ): IWolj (j2003|) and MCFOST described by 
Pinte et al. J200d . |2009[ ). Both codes are three-dimensional, radiative transfer codes using 
the Monte-Carlo method and NextGen stellar atmosphere parameters. The only significant 
differences in the modeling with each code were the range of parameters examined and a 
few of the initial assumptions described below. The goal of this modeling was to understand 
what disk parameters we are most likely to be able to constrain in the future, and what 
assumptions are most critical to define, perhaps by future observations. For example, we did 
not attempt to fit the exact shape of the silicate feature, hence the exact silicate mineralogy 
of the dust in the disk atmosphere. The integrated flux in the model feature is close to the 
data, and the thermal budget of the disk is therefore also correct and sufficiently accurate 
for this study. More detailed modeling will follow when our full BD sample is available. 

Both modeling codes parametrize the disk structure geometrically in similar ways. Typ- 
ical parameters that were fixed in our initial modeling include: the slope of the grain size 
distribution, dn{a) oc a~^'^da, minimu m grain size, and th e dust grain properties, typically 
astronomical silicates as described by iDraine fc Led (119841 ) with varying amounts of amor- 
phous carbon. In the MC3D models the surface density power law was set to oc r"^ with 
p = 1, while with MCFOST a range for the exponent from 0.1 to 1.5 was explored, depending 
on the object. The stellar parameters were fixed based on previous studies referenced below, 
but models with some variation in parameters were tested because of uncertainties in spec- 
tral types, luminosities, and stellar radii. Typical parameters that differed between models 
included: maximum grain size, total gas+dust disk mass (with an assumed gas-to-dust mass 
ratio of 100), inclination angle, inner radius, and the disk scale height parameters, Hq and 
7, h{r) = hoi^r/roy. The outer radius was typically fixed at 50-200 AU, but the SED's are 
quite insensitive to this choice, as illustrated below. 

For each object Figure [2] shows two example good-fit models whose parameters are 
listed in Table |3l We found that many of the model parameters are quite unconstrained or 
strongly dependent on other parameters or on our choice of such basic values as the exact 
stellar properties and dust properties. For example, there are the expected degeneracies 
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between stellar luminosity versus extinction, inner disk radius versus minimum grain size, 
and between disk mass and grain optical and size properties. 



4.1. SSSPM1102 



SSS PM1102-3431 was identified as an M8.5 brown dwarf in the TW Hya associa- 
tio n bv IScholz et al.l (|2005|). Photornetry a nd s pectroscopy wi t h Spi tzer were reported 



by iRiaz &: GizisI ( 120081 ): iMorrow et al.l ( 120081 ) and iLuhman et al.l (120 lOf ). and these results 
are shown in Figure 121 The object has a well-determined parallax distance of 55.2 pc 
(ITeixeira et al.l 120081 ) and no reported companions. 



The two models in Figure |2l illustrate the effect of adding 10% amorphous carbon to a 
nominal silicate grain composition, providing a better fit to the lack of a strong silicate feature 
in the Spitzer IRS data. Although most of the disk parameters are not well constrained, we 
find a fairly robust lower limit to the disk mass of a few xlO~^ Mq which is driven by the 
strong far-ir emission detected. The probability distribution is quite fiat above 10~^ Mp 
The disk inclination is also likely greater than 60°. 



0- 



4.2. IS0138 



IS0138 was identified by I Gomez fc Persil (12002) as an M 5 .5 bro wn dwarf in the Chamaeleon 
I association, though subsequent spectroscopy by LuhmanI (2004 ) f ound a spectral type of 
M6.5. Spitzer photometry of IS0138 repo rted bv iLuhman et al.l (l2008al ) and IRS spec- 
troscopy over the lO/xm silicate feature by Pascucci et al. ( 20091 ) are show n in Figure El 



Based on its presumed membership in Cha I, we assume a distance of 160 pc (ILuhman et al. 
2008bl ). and IS0138 also has no reported companions. 



Unlike the other two objects, IS0138 exhibits a weak silicate emission feature (jPascucci et al 



20091 ). The two model fits shown in Figure [2l illustrate the effect of silicate composition on 
the shape of the emission feature. The fitted inclination angle and were strongly de- 
pendent on the assumed stellar parameters. The only fitted parameter for which we have 
a reasonable constraint is the mass for which a broadly-peaked probability distribution was 
found around M~ few xlO~^ M0. The far-ir fiux implies a lower disk mass limit of a few 
xlO~^ Mq. The inner disk radius is also constrained to the likely range 0.03-0.15 AU. 
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4.3. 2M1207 



Gizid (I2OO2I) identified 2M ASS 1207334-393254 a s a likely M8 brown dwarf in TWA as 



was subsequently confirmed by IChauvin et a. 



52.4 pc has been derived by iDucourant et al 



(120041) ■ A trigonometric parallax distance of 



(120081) ■ Spitze r photometry and spectroscopy 



are described by iRiaz fc GizisI (|2008[ ) ; iMorrow et al.l (|2008[ ) and are also shown in Figure 



m 



Riaz fc Gizis (j2008l ) plot a 70 /im MIPS measurement in their Figure 6, but no fiux 
yalue is described in the text other than an uncertainty of ~ 0.4 mag, suggesting a 2a 
detection. Unlike the other BD's, 2M1207 has a clear l y iden tified companion at a projected 
separation of ~ 770 mas (55 AU) (jChauyin et al.ll2004i 120051 ). The mass and model-fit to the 



companion's SEP are a subject o f some controversy in the literature, e.g. iMohanty et al. 



( I2OO7I ) ; iMamajek &: Meyerl ( 120071 ) ; ISkemer et al.l (l201l[ ) , but for the purposes of our modeling 
we assume the primary has a substantially higher mass and luminosity than the secondary, 
and provides most of the dust heating. 

The two models shown in Figure [2] illustrate the negligible effect on the SED of trun- 
cating the outer radius of the disk due to possible effects of the known companion which 
would be expected to limit the outer disk radius to ~ 1/3 the component separation. Like 
SSSPM1102, most disk parameters are quite unconstrained by our models, but the mass 
and flaring index do exhibit peaked probability distributions. Our modeling suggests a likely 
disk mass ~ 10~^ Mq with a lower limit of a few xlO^^ Mq, though masses up to 10~^ Mq 
can also produce reasonable flts. The flaring index 7 is fairly well constrained to the range 
1.1 < 7 < 1.15. 



5. Discussion 



The strongest conclusion from our modeling is that the addition of the new Herschel data 
can provide important lower limits to the disk masses. Strong upper limits depend on as yet 
unavailable photometry at longer wavelengths and resolved imaging. With respect to most 
of the disk parameters, quite simple assumptions about the disk properties provide good flts 
to the observed SED's, but changes to different parameters produce corresponding changes 
in other fltted parameters that still leave a model with a good flt. Our derived disk masses, 
for example, suggest that the SSSPM1102 disk is the most massive, but as for all modeling, 
this conclusion depends on assumed grain properties. With such assumptions, the range 
from the least to the most massive is likely to lie between a few x 10~^ M0 up to perhaps as 
much as 10"'^ Mq (gas+dust), though with the most conservative uncertainty estimates, this 
range could be 3 times larger or smaller. The disk scale heights, Hq at 100 AU, likely lie in the 
range of 5-20 AU with modest flaring indices, 7 < 1.25, both values that are typical in disk 
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models for more massive objects. Interestingly the nominally more evolved objects may have 
less flared disks, based on their differences from IS0138, but more observations of additional 
BD's are needed to determine if this is a general effect. Es timating sub-stellar masses is 
notori ously uncertain, especially for very young objects, e.g. iBaraffe. Chabrier &: Gallardo 
( I2OO9I ). But using the spectral types mentioned above and ages of 1-2 Myr for Cham I and 
8 Myr for TWA and commo n evolutionary BP models, these three objects probably have 
masses of a few x 0.01 Mq (IBurrows et al.l 119971 : iBaraffe et al.l 120031 ). This would imply a 
ratio of disk-to-stellar mass in the range 10~^ up to a few xlO~^ for these three BD's. An 
additional conclusion from our modeling is that it is more difficult to fit the observed SED's 
with a grain size distribution appropriate for the interstellar medium than one in which there 
has been substantial growth in the maximum grain size, but even this conclusion depends 
somewhat on other assumptions about likely parameter values. Our modeling shows that the 
lack of a silicate emission feature in the SED's of 2M1207 and SSSPM1102 can be produced 
by a modest admixture of amorphous carbon to the grain composition, though this could 
also be due to a lack of small silicate grains. 

The two of our obje cts w ith complete Spitzer IRS observations have been modeled by 
Riaz fc Gizid ( 120071 . 120081 ) and lMorrow et al.l ( l2008l ). A n umber of smal l differ ences in model 
details make a direct comparison difficult. For example, iRiaz &: GizisI (120081 ) used different 
grain size distributio ns in the midplane and atmosphere of the disk to simulate dust settling. 
Morrow et al.l (120081 ) included an expanded "wall" at the disk inner edge and characterized 
their disks with accretion parameters as well as using two grain populations. Like our 
models, though, they both assumed a grai n size distribu ti on, n (a) oc a~^'^ and were both 
unabl e to constrain the o uter disk radius. iMorrow et al.l (120081 ) do not quote a total disk 



mass 



Riaz fc GizisI (120081 ) find a larger mass for the 2M1207 disk than that for SSSPM1102, 



but they did not have available the longer wavelength 160/im Herschel data. Like our results, 
both other studies find evidence for substantial grain growth and for inclinations of >60°. 
Most importantly, none of our models required either dust settling nor an expanded inner 
wall to fit the observed SED's with reduced ^ 1-4. In general we find that it is difficult to 
constrain any of the disk parameters to the degree suggested by previous modeling without 
making assumptions that are themselves rather uncertain. 

To further constrain these models the next obvious step is to extend the SED's to A ~ 
1 mm and obtain spatially resolved images with ALMA. The brighter members of our BD 
sample are likely to be resolvable with the full ALMA array. The available model parameter 
space may also be better constrained with the highest S/N spectroscopy to best define the 
sub-stellar photospheres, in particular for IS0138 which has discrepant spectral types in the 
literature. 



- 8 - 



One other BD has bee n observed oyer an even greater span of the electromagnetic spec- 
trum, 2M04442713+2512 JBouv et al l boosh from the visible to 3.5mm. They inferred a 
fairly massive disk, ~ 10-3 Mq, though even with their complete and high S/N observa- 
tions, many disk parameters were still poorly constrained. Disk ma sses for BD 's have also 
been derived for objects with millime ter wavelength photometry by iKlein et al.l ( l2003l ) and 
Scholz. Jayawardhana fc Wood! (120061 ). They estimated masses in the range of ~ 10-^-10-^ 



Mq with differing assumptions, but clearly in the same range as studies based on more 
complete SED's. It is also ir iteresting to compare these four BD's, i.e. our three plus 
2M04442713, to T Tauri stars. [Andrews k Williamsl fcoosl . boOTh have surveyed T Tau stars 
at millimeter wavelengths where the dust is optically thin and the most accurate masses can 
be determined. They find dust masses in the range of 10"^ up to perhaps as high as 10^'^ 
Mq, depending on assumed grain sizes and properties, implying total disk masses of 10"'^ 
- 10-1 Mq with our canonical gas-to-dust ratio. So the low end of the T Tau disk mass 
distribution overlaps with the high end, 2M0444 2713, of the BD distribution w ith the very 
small sample available so far for BD's. Recently iLee. Williams fc Ciezal ( 1201 ll ) have found 
substantially lower disk masses around T Tauri stars in the 2 Myr old IC 348 association 
than in the very youngest, nearby star- forming clouds, Taurus and Ophiuchus. Since TWA 
and probably Cham I are both older than Taurus and Ophiuchus, it may be best to compare 
our BD disks to T Tauri disks in older star-forming clouds. With the completion of our data 
set within the coming year, we should have enough statistics to discern any such trends with 
better significance. 



6. Summary 

We have detected all three of our first program objects at 70/im and one at 160/im. 
These observations represent by far the most sensitive far-infrared photometry of brown 
dwarfs. Our modeling shows that the SED's can be fit with simple geometric disks that 
do not require an inner wall nor dust settling. The addition of our Herschel measurements 
provides much stronger lower limits to the masses of the circumstellar disks because most of 
the disk mass is at large radii and relatively cool. The implied disk masses are probably well 
below those surrounding many of their more massive counterparts, the T Tau stars. These 
disks are likely to be optically thin in the far-ir perpendicular to the disk plane outside 1 
AU, though optical depths through the midplane are still high, > 1000. Interestingly the 
least massive disk is found around the nominally youngest BD, IS0138 in Cham I, while the 
disks around the BD's in TWA, presumed to have an age of order 8 Myr, are probably more 
massive. Future ALMA observations will enable more accurate masses and disk sizes to be 
determined for many of the objects in our sample that are bright enough to be resolved. For 
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example, the most massive disk in this small sub-sample, that around SSSPM1102, would 
likely have a flux of slighly over 1 mJy at SBO/im, and these three objects are some of the 
fainter members of our total sample. 
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Tabic 1: Observational Parameters 



Parameter 



Value 



Comments 



AOR Type 
Wavelengths 
Number of Scan Legs 
Scan Length 
Cross Scan Step 
Scan Angles 
Repetitions 

Peak Intg Time Per Pixel 



PACS Mini-Scan-Map Two Crossed AOR's 
70/xm, 160/xm 
8 
3' 
4" 

70°, 110° Relative to Detector 

7 Per AOR 

504 sec Per AOR 



Table 2: Observations Summary (Program ID = GTl_pharve01_2) 



Object 


RA/Dcc Center (2000) 


AOR's 


Obs. Date 


70/im (mJy) 


160/im (mJy) 


SSSPM1102 


11 02 09.8 -34 30 36 


1342221849/50 


29 May 2011 


7.3 ± 1.0 


7.1 ± 1.5 


IS0138 


11 08 19.0 -77 30 41 


1342218699/700 


16 Apr 2011 


3.7 ± 0.6 


< 15 


2M1207 


12 07 33.4 -39 32 54 


1342202557/58 


10 Aug 2010 


7.0 ± 0.8 


< 7 
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Table 3. Parameters (and Range Explored) For Models in Figure [2] 





IS013S 


2M1207 


SSPM1102 


Parameter 




^ VrllllP fRrlTlffpl 


^ VrllllP fRririffpl 


stellar 








Teff (K) 


2900 


2600 


2600 


Rstar (R-o) 


0.35 (0.2-0.5) 


0.24 (0.2-0.3) 


^0.35, ^0.27 (0.2-0.4) 


Luminosity (L0) 


0.0168 


0.0046 


0.0059 


Ay 


0.5 


0.0 


'^O.O, ^1.5 


Dist. (pc) 


160 


53 


56 


Disk 








Incl. (deg) 


<60 (0-90) 


70, ^ 78 (0-90) 


^^80, ^^66 (0-90) 




0.08 (0.08-40) 


0.015 (.005-.03) 


^0.015, ^0.006 (.0035-.015) 


Router (AU) 


100 


^ 75, ^ 20 


75 


i/o(AU)@100AU 


20 (5-30) 


9 (5-15) 


5 (5-15) 


7 


1.25 (1.0-1.25) 


1.125 (1.00-1.125) 


n.05, ^1.07 (1.0-1.125) 


-P 


0.8 (0.5-1.5) 


1.0 (0.5-1.5) 


^0.5, '^0.1 (0.1-1.5) 


Dust 








2 Log{Maisk) (Mq) 


-5.2 (-6.7- -3.7) 


^-5.0 '^-5.2 (-6.5- -3.5) 


-4.0 (-6.5- -3.5) 


amm(/im) 


0.05 


0.1 


0.05 




1000 (10-1000) 


1000 (10-1000) 


10 (10-1000) 


Power law 


-3.5 


-3.5 


-3.5 


Silicate 


^0.95 DL, ^0.95 Olivine 


0.95 DL 


^^0.90 DL, ^1.00 DL 


Amorphous Carbon 


0.05 


0.05 


'^O.IO, ^0.00 



^Wliere two values are given, the first is for the "a" model and the second for the "b" model. Otherwise 
parameters were the same for both models. 

^Assuming gas-to-dust mass ratio of 100. 
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Fig. 1.— Top, left to right - 70/im images, 90" square: SSSPM1102, IS0138, and 
2MASS1207; bottom - 160/im images. Circles are 30" diameter centered on the nominal 
source positions. 
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Fig. 2. — SED's of the three observed BD's (open circles). Spitzer IRS data are shown in 
green. Sohd hnes show SED's from the example fitted models listed in Table [31 Dashed lines 
show the bare sub-stellar photospheres used in the model, extincted by the chosen A^. 
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